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ABSTRACT 

In the present work, we extend the classic halo model for the large-scale matter distri- 
bution including a triaxial model for the halo profiles and their alignments. In partic- 
ular, we derive general expressions for the halo-matter cross correlation function. In 
addition, by numerical integration, we obtain instances of the cross-correlation func- 
tion depending on the directions given by halo shape axes. These functions are called 
anisotropic cross-correlations. With the aim of comparing our theoretical results with 
the simulations, we compute averaged anisotropic correlations in cones with their sym- 
metry axis along each shape direction of the centre halo. From these comparisons we 
characterise and quantify the alignment of dark matter haloes on the ACDM context 
by means of the presented anisotropic halo model. As our model requires multidimen- 
sional integral computation we implement a Monte Carlo method on GPU hardware 
which allows us to increase the precision of the results whereas it improves the per- 
formance of the computation. 

Key words: galaxies: groups: general, dark matter, large-scale structure of universe. 



1 INTRODUCTION 



In the mid-twentieth century, man fc Scottl (| 19521 ) pub- 
lished a work which would be the basis of the Halo Model. In 
that paper the authors proposed a novel model for the spa- 
tial distribution of galaxies in the Universe. To build up their 
model they imposed four main assumptions: (i) galaxies in 
the Universe lie inside galaxy clusters; (ii) the number of 
galaxies belonging to a cluster obeys a probabilistic law; (iii) 
there is a probabilistic law which describes the distribution 
of galaxies inside such galaxy clusters; and (iv) the cluster 
centres are placed in the space under a quasi-uniform proba- 
bilistic law. Although this formalism was originally thought 
to describe the distribution of galaxies it was later gener- 
alized to describe the dark matter distribution (e.g. , ISeliakl 
l200d : iPeacock fc Smithll200d : ICoorav fc Shethll2002t ) . 

In its standard form the model assumes that all the 
matter in the Universe reside in matter haloes which 
have a spherical density profile. In addition the halo-halo 
self-correlation at large scales is governed by the linear 
power spectrum. In this sense, the standard halo model 
does not include any assumption about alignments in the 
halo-halo distribution. However, it is well known from nu- 
merical simulations, that haloes exhibit mildly aspheri- 
cal shapes, with a slight p r eferen ce towards prolate forms 
(see for instance |Paz et ail (|2006h and references therein). 
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Moreover, it has been shown that their orientations are 
related to the surro unding structures such as filaments 
and large-scale w alls (Bond ct al. 1996; Colbcrg ct al. 2005 
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iBett et al.ll2007; Kasun fc Evrarc 



iBasilakos et alj|200d; lHahn et al 



2005.: lAllgood et al 
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iNoh fc Cohn|[201ll ) . This results are also supported by 



observations showing a good agreement with those obtained 
from the g alaxy group shapes alignment with the galaxy dis- 
tribution (iPaz et alJ|200S lLau et al J [20111 ; |Paz et alJ^Olll ; 

"Hi 



ISmargon et al. 20121 ) . Most of these studies perform statis- 
tics on the inclinations of halo axes with respect to directions 
defined by the surrounding structure. This structure is char- 
acterized by the particle distribution through various topo- 
logical sigjiatures_^uchasfilaments, walls or voids (see for in- 
stance [CauturTeTlil] [2013]). While these studies agree about 
the presence of alignments within the ACDM context, the 
magnitude of the alignment effect depends strongly on the 
geometrical definition of the surrounding structure. On the 
other hand. lPaz et al.l (|2008l ) have introduce a new approach 
which is independent of any definition of surrounding struc- 
ture. This approach uses a modified version of the classical 
two-point cross-correlation function to quantify the align- 
ment between centre objects and the Large Scale Structure 
(LSS). Furthermore, this modified cross-correlation func- 
tion allows a robust comparison between numerical sim- 
ulations and observational data. This approach has suc- 
cessfully characterized alignments between galaxy shapes, 



2 Sgro, Paz & Merchdn 



group shapes and simulated dark matter haloes with sur- 
roundi ng structure traced by galaxies or simul a ted struc- 
tures llPaz et al.ll2008l: iFaltenbacher et a.1.1 120091: IPaz et all 
l201ll ; IFaltenbacher et alj|2012i : ISchneider et al.||2012fl. 



lies inside haloes, the density field p(x) can be computed as 



Relevant to this work are the studies by Smith fc Watts! 



(|2005l ). These authors investigated the effects of halo triax- 
iality and LSS alignment on the isotropic power spectrum. 
They found that the effect of halo triaxiality manifests a 
small suppression of power at the 1-halo term, whereas in the 
2-halo regime they found an upper limit of 10% of clustering 
increment, by imposing perfectly aligned haloes. These 
correlations on halo triaxiality can affect cosmic shear stud - 
ies used to analyse weak lensing ()Bridle fc Abdallal l2007f l . 
Moreover, an intrinsic alignment in the structure traced 
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In this paper we derive an analytical model for the 
anisotropic halo-matter cross correlation function. We im- 
plement an extension of the classical halo model for the 
large-scale matter distribution, that includes tria xial mod- 
elling of halo profiles (following I Jing &~S uto 2002) and their 
alignments (described in ij2]) . The use of the correlation func- 
tion instead of the power spectrum allows a straightforward 
implementation of alignment effects and triaxiality of the 
haloes. Furthermore, this anisotropic cross-correlation func- 
tions can be easily measured from observat ions and numer- 
ical r esults as showed in previous works |Paz et al.l 120081 . 
l201ll ). Therefore, the main goal of the present work is to 
obtain a suitable analytic model for these functions in order 
to compare with numerical. 

The organization of this paper is as follow. In Sj2]we in- 
troduce the anisotropic halo model and we describe some 
assumptions performed. In 331 we describe the numerical 
implementation. In f|4] we compare the halo-matter cross- 
correlation functions obtained from our model with those 
estimated from N-body simulations. Finally, we discuss and 
summarize our results in SJU 



2 ANISOTROPIC HALO MODEL 

In this section we introduce the analytic development of 
a new halo model which takes into account the triaxial 
nature of the haloes and its alignment with the environ- 
ment. As in the standard halo model, we will assume that 
the matter distribution is made up of distinct haloes in 
a wide range of masses. We characterize each dark mat- 
ter halo in terms of its mass normalized density profile 
U(f,m,a,e) = p(r,m,a,e)/m where r is the position re- 
spect to the centre of an halo of mass m, the components of 
a axe the eigenvalues of the shape tensor of the halo (with 
0.1 > 0,2 > 0,3), and <f indicates the direction of the eigen- 
vectors (i.e. its components are the tree Euler angles of the 
shape ellipsoid orientation). Assuming that all dark matter 



p(x) = y~] rriiU(x — Xi,mi,ai,6i) 



(1) 



where N is the total number of haloes in the Universe. Sim- 
ilarly, the number density of haloes n c (x) can be expressed 
as follow 



z{x) = S(S - Xi)il>(mi, Si, 



(2) 



where the factor i/j(m,a, e) is the selection function, and 
8(x — Xi) is the Dirac delta function centred at the halo 
position. 

The probability distribution p(m, a, e) of haloes with 
mass m, shape vector a and orientation e can be computed 
using the equation 

IN \ 

p(m, a, e) = / ^ 8(m — rrii)5(a — Si)5(e — ?i)8(x — Xi) \ , 

(3) 

where the angle brackets denote the average over an ensem- 
ble of realizations. In case that shapes and orientations are 
not taken into account, this function reduces to the stan- 
dard mass function. From the above definitions, the mean 
dark matter density and the mean number density can be 
computed as: 



dm dade mp(m, a, e) 



n c = I dmdddep{m,a,e)^l){m,d,t) 



(4) 



(5) 



Hereafter the integration limits should be understood as 
taken over the whole range where the integrand is valid. The 
^(m, a, e) function allows us to impose some restrictions on 
the mass range, shape and orientation of the centre haloes. 

Using these function we are able to compute the halo- 
matter cross correlation function as follow: 

icQg) _ A/ p{x + r) _ ^ 
— c J \ P 
' n c (x)p(x + r) n c (x) p(x + r) 



n c p n 
n c (x)p(x + r) 



+ 1 



71c p 



(6) 



where we have assumed the ensemble average can be re- 
placed by an space average: 



and similarly 



Replacing [T] and [2] in equation [6] we obtain: 
£,h m {f) = -1 + =!= ( (y^5(x- x^tptrrir, Oi, u) ] 

prtc \Wi / 



x \^mjU(x + r- Xj^j^j,^)^ (9) 
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It is possible to split this equation in two terms, those 
where i — j (hereafter 1-halo term), and those where i 7^ j 
(hereafter 2-halo term). Consequently we define: 

1 / N 

x miU(x + f— Xi,rrii,ai,ei) ) (10) 



dKf) = - 1 + 



l / fJL 

nc We? 

y m j U ( x + f - x j • , m , • , a j , e ) ■ ) j 



(11) 



Rewritting equation [10] as: 
11^ 

ffcm(r) = ^77 y~]rniifi(mi,ai,ei 
pn c V 



x y dxS(x — Xi)U(x + f—Si,mi,ai,ei) 



1 1 

77 ' rnixj){mi,ai,€i)U{r,mi,ai,ti) 

pn c V f-f 



dm dade mip(m, a, e)U(r, m, a, e) 



1 ^ 

— ^ 5(m - rrn)S(a - ai)<5(e — u 



V 

Using equation [3] the final expression for the 1-halo 
term results: 



(12) 



£i,„Jr) = / dmdad€mU(r,m,a,e) 

pn c J 

x p(m, a, e)tp(m, a, e). 



A similar procedure can be used to derive the 2-halo 
term from equation 1111 

t^hinir) = — 1 + / dmi dd\ del dm 2 dd 2 de* 2 dy 

pn c J 

x m 2 U(f- y, m 2 , 3 2 , e 2 )^(mi, 3i,ei) 
^ N N 

x — ^2 3( mi ~ •m i )8{di - Si)5(ei - ei) 



j=0 i=l 



x <5(m 2 - mj)5(a 2 - dj)S(e 2 — £j)5(y - (fj - a?j)) 

(13) 

The joint probability to have a pair of haloes with a 
given set of properties (mass, shape and orientation) at a 
distance \y\, can be written 

AT AT 

i=o i=i 

x 5(m 2 — rrij)S(a2 — 3j)5(e*2 — ej) 

x 8(y- (xj - Si)) = pip 2 (£i, 2 (y) + 1), (14) 

where — p(mi,3i,€i), i = 1,2, are the probabilities to 
have a halo with a given set of properties (mass, shape and 
orientation), and £1,2(1/) + 1 = £(2/1,2, mi, Si, el, m 2 , a 2 , e 2 ) + 



1 is the join probability to have a pair of haloes with the 
given properties separated by a distance 3/1,3. In the right 
hand side we have assumed that the probability of halo prop- 
erties (pi and p 2 ) and the join probability to have a pair of 
haloes separated by a distance y are independent. Therefore, 
using this last equation, the final expression for the 2-halo 
term is given by: 

£hm(r) = / dmi dai del dm 2 da 2 de* 2 dy 

pri c J 

x miU(r- y, m 2 , a 2 , t2)piP2^i,2(y) 

x tp(mi, 3i, el) (15) 

Unlike standard models, equations[l2]and[l5]give a halo 
model for the two-point anisotropic cross-correlation func- 
tion taken into account the shape and the orientation of 
haloes. We will devote the next two sections to explain our 
implementation of this model in order to compare with the 
corresponding measured correlations in the simulation. 



2.1 The 1-halo term 

In order to estimate the 1-halo term (equation [12} , we need 
to adopt models for the halo normalized density profile (U) 
and for the distribution of halo parameters (p). Regarding 
the density profile , we use the triaxial model described by 
IJing fc Sutof l|2002l ) (hereafter JS). This model assumes that 
isodensity surfaces of dark matter can be described by tri- 
axial ellipsoids with a shape vector a and an orientation e*. 
A simple way to parametrize tho se surfaces is through the 
radial parameter R as defined bv lSmith fc Watts! (|2005l ): 



«3 2 



2 2 2 

z y x 

9 "T" 9 ' '1 

Oui z Ou2 Z £13 



(16) 



Accordingly to JS, the former definition all ows to express 
the h alo density profile similarly to that in iNavarro et al.l 
(|l997h (NFW profile): 



5c 



p(R) 

Pent (R/Ro)(l + R/Ro) 



(17) 



where S c is a characteristic density and Rq is a scale radius. 
Analogously to the scale radius in NFW profile, JS define 
Ro = c e /r e , where c c is the concentration parameter and 
r e is the characteristic radius. JS also provides the relations 
c e — x c vir, T e = X r vii', where Cvir, r vlr are the virial concen- 
tration and radius, respectively. It is possible to relate the 
characteristic density S c with the concentration parameter 
by requiring that the mean density inside an ellipsoid with 
major axis r e to be A e times the critical density, obtaining: 



5 C = 



3 ln(l + C e ) -Ce/(1+Ce) 



(18) 



JS provides a relation involving this the equivalent density 
A e as function of the virialization density in a spherical 
NFW: 



5A V 



(l Ml 2 



(19) 



A value of x = 0.45 was obtained by the authors fitting 
numerical simulations. We found that a \ value of 0.5, is 
also consistent with simulation data and provides a better 
agreement between the NFW profile and the JS profile at 
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the spherical shape limit. In this limit case, our choice for 
the x parameter gives similar values for the 8 C and the char- 
acteristic density of the NFW profile in a mass range of 10 8 
to 10 15 M Q . 

Once the normalized density profile has been estab- 
lished, we must set up the probability distribution p(a, m, e). 
For simplicity, we assume that the orientation e*and the mass 
m of a given halo are independent random variables, whereas 
the halo shape vector a is only dependent on the mass. Con- 
sequently the probability distribution can be written as: 



p(m,e,a) = n(rn)p(e)p(a\m) 



(20) 



where n(m) is the halo mass function, p(e) and p(a\m) are 
the orientation and shape probability distributions, respec- 
tively. For instance, in case of a uniform probability for the 
halo orientation on the sphere, the density function p(e) 
takes the form: 



p(e)de = p{a, j3,^)dadj3d"j 
111 

= dad(cosB)d r y 

27T 2 27T K Hl ' 



(21) 



where the Euler angles are restricted to the ranges: 
> a ^ 2tt, ^ /3 ^ 7T, and > 7 > 2?r. 

Given that equation [16] can be written in term of the 
ratios 021 = 02/01 and 032 = 03/02, we define the p(a\m) 
conditional probability distribution as a function of these 
quotients p(a2i, 032 where both arguments are restricted 
to the interval (0, 1). With this simplification, we are able 
to exchange da — > da2idaz2 in equations [12] and 1151 We 
approximate this distribution by a product of two Gaussian 
laws and quadratic factors that ensure null probability in 
the interval limits 0. Consequently, the shape distribution 
takes the following form: 

(q 21 -4) 2 

p(«21, tt32|m) OCd21 (1 — <X2l) e 2^ 



(°32-B) Z 

xa32(l-ci32)e 2^2 



(22) 



where A and B are two mass dependent parameters which 
must be estimated, whereas a fixed value a = 0.1 is 
ad opted. A theoret ical approach to this function was made 
bv lLee et all (|2005h . However, their results are not in good 
agreement with the numerical simulation. On the other 
hand, JS used a direct measurement of the axis ratios on 
dark matter haloes extracted from N-body simulations with 
the aim to find the density distribution. In our work, we 
adopt the analytical form given by equation 1221 which re- 
sembles the shape distribution obtained from the simula- 
tions. The two free parameters of the adopted distribution 
are used to fit the cross-correlation function. 

The last ingredient required in order to estimate the 
1-halo term is the mass function n(m), which describes the 
abundance of dark matter haloes in a mass interval around 
m. This fun ction is estimated u sing the analytical formulae 
provided bv lSheth etal] (|200ll ). 



2.2 The 2-halo term 

In the current section we will set the corresponding func- 
tions for the 2-halo term. To compute equation [T5] we need 
to specify the function £1,2(2/) = £(j7i,2|?7ii, 0*1 , <fi, 7712, 0*2, £2). 
On large scales the relation between haloes and matter over- 
density is deterministic, therefore, as in the standard halo 
model, we can approximate £1,2 by the following equation: 

f;(y)i,2 ~ /(yi,2|ai,ei,02,€* 2 )6(w.i)6(m2)6m(yi,2) (23) 

where fe(mi) and 6(7712) are the bias factors for the centre 
and tracer haloes, respectively, and £u n {yi,2) is the Fourier 
Transform of the Linear Power Spectrum (LPS). We have ex- 
plicitly introduced the factor / (2/1,2 \Si , ei , 0*2, £2) to take into 
account the different alignments. We have adop ted the ana- 
lytica l expression of the bias function given bv lSheth et al.l 
(2001) whereas the LPS i s com puted using the approach 
given bv lEisenstein fc Hul (|l999h . 

The / function defined above can be used to describe 
two types of alignments: (i) the alignment between the po- 
sition of tracer haloes and the major-axis direction of cen- 
tre halo, (ii) the alignment between the orientations of cen- 
tre and tracer haloes. It is reasonable to expect that the 
last kind of alignment has a second order effect on the 
model. This can be explained by t he fact that our model , 
and previous works (see for instanc e van Daalen et aT]|2012l : 
ISmith et aD 120061 : ISmith fc Watts! I2005T ). are insensitive to 
the halo shapes in the 2-halo term. This is expected given 
that the difference between the major and minor axes of a 
halo is small compared to the distances, at these scales, be- 
tween centre and tracer haloes. Consequently, for the sake 
of simplicity, we have assumed that the alignment function 
does not depend on the alignment between the major axes 
of centre and tracer haloes. Furthermore, we will consider 
this function as independent of halo shapes. Based on the 
previous considerations we adopt the following expression: 



/(yi,2|ei) oc e Tc? + De 



(24) 



where 9 and <j> are the polar and azimuth angle of 7/1,2 in a 
spherical coordinate system defined by the directions given 
by ei, respectively. The first (second) term in the preceding 
equation can be understood as the excess of correlation along 
the major (intermediate) axis. As the alignment function is a 
probability density distribution, the constant of proportion- 
ality is defined by its normalization. The C and D factors 
are free parameters that we will determine in EI4.2I Figure 
[T] shows an illustrative example of the alignment function 
with parameters C = 0.3 and D — 1.0. It should be noticed 
that even when 6 approach to and <j> tend to be undeter- 
mined, the alignment nicely becomes constant by construc- 
tion, avoiding any undesirable behaviour on its evaluation. 
The increment of the function towards low 9 values repre- 
sents the alignment of the matter distribution around the 
major shape axis, whereas the increment towards 9 = n/2 
and = represents the corresponding alignment with the 
intermediate axis. 



3 IMPLEMENTATION 



1 These factors ensure that a perfect sphere (021 = 032 = 1) or 
an ellipsoid with a null semi-axis have a null measurement 



In this section we first describe the numerical method used 
to compute the integrals ll2l and [T5l Secondly, we show the re- 
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Figure 2. Results of the anisotropic cross-correlation functions along the halo shape axis directions. Centre halo masses span a range 
from 10 12 to 1O 16 M0/i _1 . Solid, dotted and dashed lines correspond to major, intermediate and minor axis direction, respectively. Each 
panel corresponds to a different set of parameters as explained in the text. 



following expression: 



/ = 



f(x)dxid,X2---dx 71 



(25) 




Figure 1. Normalized alignment density distribution around 
dark matter haloes with parameters C = 0.3 and D = 1.0 (see 
equation 1241) . 



suits for the three-dimensional anisotropic cross-correlation 
function. 



3.1 Integration method 

Given the high order integrations involved in the 1-halo and 
2-halo terms of our model, standard numerical methods (e.g. 
quadrature rules) are not feasible to implement. In order to 
compute the multidimensional integrals we have employed 
a Monte Carlo technique (MC). According to this method, 
the m— dimensional integral 7 of a given function f(x) in 
the volume V C R m can be approximated by means of the 



where Xi = (xi , .., a;™), are TV uniform random vectors 
within the computational volume V c D V. It can be shown 
that the uncertainty in the determination of I decreases as 
7V~ 1//2 independently of the dimension of the integrand. It 
is possible to improve the precision of the MC integration 
by using the importance sampling method. This method 
requires the specification of a density function w(x) for the 
random variable x. Thus, in the right hand side of equation 
[231 f{xi) must be replaced by f(xi)/w(xi). It is worth to 
mention that using standard integration methods, for a 
fixed number of subdivisions N of the volume V, the error 
erj quickly increases as the number of dimensions increase 
(pi cx N~ 1/m ). 

As stated in the preceding paragraph, to estimate the 
integrals [12] and [15] we should generate as many random 
numbers as possible. The Monte Carlo integration is an 
embarrassingly parallel problem, so we have employed the 
Nvidia CUDA extension for C language in order to compute 
the integrals. This choice allows us to increase the amount of 
random numbers and consequently to enhance the precision 
of the results while the total running time grows moderately. 



3.2 The anisotropic tridimensional 
cross-correlation function 

Before comparing our halo model with results obtained from 
numerical simulations, we compute the cross-correlation 
function along the shape axis directions of the centre halo 
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in order to give a qualitatively analysis of the parameter 
dependence. For this purpose, we fix the orientation of the 
centre halo such that the directions of the shape axes corre- 
spond to cartesian axes (denoted by ei, ea, respectively). 
Accordingly to this choice, the integral in equation [12] l)15p 
is reduced to a 3-D (12-D) integration. We define three cor- 
relation functions by taking r || e* where i = 1, 2, 3. 

With the aim of describing qualitatively the behaviour 
of our model for different values of A, B, C & D, figure 
[2] shows the anisotropic cross-correlations for haloes in the 
mass range 10 12 — 10 16 Mq/i" 1 . The three directions, major, 
intermediate and minor axis are shown with solid, dotted 
and dashed lines, respectively. Each panel corresponds to a 
different set of parameters. Upper panels show the correla- 
tion function behaviour when the shape parameters A and 
B are changed. As can be seen in panel (I) when A = B and 
C & D are properly set, the three functions show roughly 
constant differences. Panel II (III) shows the results for pro- 
late (oblate) halo shapes B < A (A < B) whereas panel 
IV shows the behaviour for a nearly spherical halo shape, 
i.e. A = B w 1. Regarding to the alignment function, lower 
panels in figure [2] illustrate the effect of changing C and 
D parameters. Panels V and VI (VII and VIII) show the 
influence when assuming low and high values of D (C), re- 
spectively. As can be appreciated, for large values of the 
C parameter the alignment vanishes because the alignment 
distribution function is nearly constant. From this qualita- 
tive analysis, it can be seen that the behaviour of the model 
shows a large versatility with the parameters variation. Fur- 
thermore, as expected, the parameters A & B mainly affect 
the 1-halo term whereas the parameters C & D are involved 
only in the 2-halo computation. 



4 COMPARISON WITH NUMERICAL 
SIMULATIONS 

4.1 Cosmological Simulation and estimation of 
the properties of the dark matter haloes 

In order to test our model, we have used a collisionless 
simulation of 1024 3 particles covering a periodic volume 
of 500 3 (h~ 1 Mpc) 3 . The initial conditions at redshift 
~ 50 were calculated assuming a spatially flat low-density 
Universe with cosmological parameters taken from the 
Seven- Year WMAP results (|Larson et alJ 1201 ll ) (matter 
density fi m = 1 - fi A = 0.272, Hubble constant H a = 70.2 
km s _1 Mpc -1 , and normalization parameter as = 0.807). 
With these parameters the resulting particle resolution is 
m p = 8.78819 x 10 9 h' 1 M Q . The run has been performed 
usi ng the second version of the GADGET code developed 
bv lSpringell (|2005l ). 

The identification of particle clumps was carried out 
by means of a standard friends-of-friends algorithm with a 
percolation length of I = 0.17 £'~ 1//3 , where v is the mean 
number density of DM particles. For this study we only 
kept dark matter haloes with at least 20 particles. 

As usual, halo shapes a and orientation e are defined by 
means of eigenvalues and eigenvectors of the shape tensor 
obtained from halo particle distributions (see for instance 



iPaz et all l200d l201ll ). The anisotropic cross-correlation 
function is then defined relative to this h alo shape axes, b y 
following the same procedure described in I Paz et all (|201lf ). 
This method computes cross-correlation functions by count- 
ing halo-particle pairs in conical volumes around the eigen- 
vector directions. The angles of these volumes are selected 
such that they are mutually disjoint. In order to compare 
these correlation functions with our analytical model, we 
compute the average anisotropic halo model over equiva- 
lent volumes. To this end, we have computed the correlation 
functions by averaging our halo model over a conical section 
around the three shape axes directions. This average is per- 
formed over the two angular components of r on equations 
1121 and [15] increasing the number of integrals by a number 
of two in each equation. 

At this point, we are ready to compute the anisotropic 
halo model for a given set of A, B, C and D parameters. In 
the following subsection we estimate this set of parameters 
through a maximum likelihood method. However, we will 
first try to recover the isotropic classic halo model by aver- 
aging over a spherical volume (obtained by taking a conical 
volume with angle equal to n). Figure [3] shows the results 
for this computation as a function of centre-halo mass rang- 



ing from 10 



10 Mq on logarithmic intervals of 0.5. 



The 1- and 2-halo integrals are showed separately. The mea- 
surements of the cross correlation function on the numerical 
simulation described above are showed as filled circles for 
the same mass ranges. Jacknife errors are represented by 
the shaded grey bands. On the left hand panel, the models 
are obtained setting A = B = 0.99 and C = D = 1.00. With 
this choice for A and B parameters, dark matter haloes have 
roughly spherical shapes. The right side panel shows the re- 
sults when the parameters A and B are set to their best 
value, estimated using the minimization method described 
below. As can be appreciated the differences between the 
simulation and both models are more important at the scales 
corresponding to the 1-halo term. Additionally, the model 
with triaxial haloes describes better the numerical results 
than the model with spherical haloes. 

We have found that the model with triaxial haloes has 
an average discrepancy of %5 with the results obtained from 
numerical simulations, whereas the model with spherical 
haloes has an average discrepancy of %20. It should be no- 
ticed that the discrepancy between the simulation and the 
models is more significant for high mass haloes. This is an ex- 
pected behaviour since the haloes have more non-spherical 
pr ofiles. This is in quali tative ly agreement with the w ork 
of ISmith fc Watts! d2005h and Ivan Daalen et all (|2012T ). In 
the first paper jSmith fc Wattsll2005r ) the authors computed 
the classic power spectrum of matter, finding a small but 
noticeable effect of halo shapes at large frequencies (small 
scales). On the other ha nd, the authors of the second work 
( van Daalen et all I2012T ) have found an effect of %20 on 
the galaxy correlation function at small scales when satel- 
lite galaxies are sp hericalized ar o und t he central galaxy. In 
the same direction, ISmith et all (|2006l ) found a similar be- 
haviour on the bispectrum. They show that the use of triax- 
ial haloes in the bispectrum model results in a suppression of 
~ 7% on scales k > O^hMpc" 1 relative to a model based on 
spherical haloes. Moreover, as our anisotropic model, their 
approach is insensitive to halo shapes on large scales. 
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Figure 3. Isotropic Halo Model. Filled dots in both panels show the cross correlation functions measured on the numerical simulation 
for different mass ranges expanded from 1O 12 M0 to 10 15 Mq. Solid lines show the results obtained by computing the anisotropic halo 
model averaging over a conical volume with angle equal to n. The 1-halo and 2-halo terms are showed separately. On the left panel the 
model results are obtained by setting A = B = 0.99 and C = D = 1.00, whereas, in the right panel, A and B parameters correspond 
to the best values obtained following the \ 2 Minimization Method. As can be appreciated the model obtained using the best values 
describes better the numerical results than the model with spherical halo profiles. 



4.2 Estimation of parameters: Maximun 
Likelihood Method 

Before computing and comparing the anisotropic halo model 
with the simulations, we need to determine the A, B, C 
and D parameters. To this purpose we use a standard \ 2 
minimization method. Given that the parameters A and B 
affect principally the 1-halo term while C and D are only 
involved in the 2-halo term calculation, we explore the pa- 
rameter space with two sets of two parameters, separately. 
We determine these parameters by fitting, on the one side, 
the 1-halo term {A and B) and, on the other side, the 2-halo 
term (C and D). This is performed by fitting each halo term 
over two disjoint scale intervals separated to avoid mixture 
between the both terms (see light grey vertical stride on 
figure [5} . 

To estimate the best fitting parameters we apply the 
Maximum Likelihood Estimation method. Figure 0] shows 
the likelihood levels of the parameter estimation for three 
different mass ranges, 10 12 -10 12 ' 5 , 10 l:! ' 5 -10 14 and 10 14 ' 5 - 
10 15 , from left to right respectively. Upper panels show the 
results for A and B parameters while lower panels show the 
likelihood levels for C and D parameters. The shaded iso- 
contour regions in black, grey and light-grey colour, show 
the 1 — <t, 2 — a and 3 — a levels for the likelihood function, 
respectively. The filled white circles on the black region in- 
dicates the best fitting set of parameters. 

As expected, most massive haloes have a less spheri- 
cal shape and more prolate mass profile. This can be seen 
on the upper panels of the figure [4] where more massive 
haloes show lower A and B parameters, meanwhile the ratio 



A/B fairly increases. These results are consistent with cur- 
rent agreements about halo shapes and structure formation 
II Paz et al.l 120061 : lAltav et al ] 120061 : iRobotham et al 1 l-'IHI.4 
lLau et alj|201lf ). It is worth to mention that the uncertainty 
in the parameter space growth with the halo mass. This can 
be understood in terms of the goodness of our model to 
fit the simulations. The adequacy of our model at the 1- 
halo regime is better for low mass haloes than for high mass 
haloes. 

The lower panels of figure [4] show the behaviour of the 
C and D parameters estimation. It can be seen that while 
C decreases with the halo mass, the D parameter increases. 
On the other hand, the 1 — a region is smaller for high than 
for low mass haloes. The shape of the isocontours on the 
left panel is elongated along the D parameter. This is due 
to the weak dependence of the alignment function on the D 
parameter for large C values. 

4.3 Anisotropic halo model applied to the 
simulation 

Using the best parameter values estimated above we com- 
pute the anisotropic halo model for the same three mass 
ranges defined in the previous subsection. The results are 
showed on figure [5] where each panel corresponds to low, 
intermediate and high mass sample from top to bottom. 
Dotted lines in the upper part of each panel show the 
anisotropic cross-correlation functions estimated from the 
numerical simulation. Each of these measurements are com- 
puted adopting a semi-angle equal to 7t/4 for the conical vol- 
ume. Correlation function measurements along the major, 
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Figure 4. Likelihood levels of the parameter estimation for three different mass ranges, 10 12 — 10 12 ' 5 , 10 13 ' 5 — 10 14 and 10 14 ' 5 — 10 15 , 
from left to right. Upper panels show the results for A and B parameters while lower panels show the likelihood levels for C and D 
parameters. 



intermediate and minor shape axes are showed in squares, 
circles and diamonds, respectively. Dark grey shaded bands 
around these values correspond to the Jacknife errors at 1 — a 
level. Anisotropic Halo Model results are displayed in solid, 
dashed, and dot-dashed lines, corresponding to major, inter- 
mediate and minor axes, respectively. At the bottom of each 
panel, we show the quotients of each correlation function 
along shape axis directions to their corresponding isotropic 
cross correlation function for the model and the simulation. 
The light grey vertical bands indicate the scale ranges used 
to estimate the best fitting parameters. As can be seen, for 
the three alignment functions, there is a good agreement 
between anisotropic halo model predictions and simulation 
measurements. 

It is worth to mention that the predictions of our model 
describe better the simulation results in the 2-halo than in 
the 1-halo regime. This could be arise from the fact that sim- 
ulated dark matter profiles become more spherical towards 
the centre, whereas our adopted JS profile assumes constant 
shape. This can be appreciated in lower box of each panel 



in figure [5] where the quotients of simulation measurements 
become smaller at small scales while for our models remain 
roughly constant. On the other hand the predictions of our 
model describe successfully the simulation results on the 2- 
halo regime over a broad scale range. 



5 CONCLUSIONS 

In this work we have developed a formalism to compute the 
anisotropic halo-matter cross-correlation function. Main in- 
gredients added to the classic halo model are the triaxial 
nature of halo profiles, their probability distribution on axis 
ratios and the probability of alignments between two haloes. 
Since we are interested to quantify the anisotropic alignment 
of dark matter haloes through the cross-correlation func- 
tions, we have chosen to develop this formalism on real space 
opposite to the standard halo model which is computed on 
Fourier space. This choice impose the calculation of multidi- 
mensional integrals in both, 1-halo and 2-halo terms, which 
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Figure 5. Anisotropic Halo Model for centre haloes in the mass 
range 10 1200 - 10 125 , 10 13 - 50 - 10 14 and 10 14 so - 10 15 ' from 
up to down, respectively. Solid lines show the results obtained 
by our model with parameters set to the best value estimated as 
described in the text. The scale ranges used for these estimations 
are indicated by vertical light grey strides. Black points show the 
measurements of the cross-correlation from numerical simulation 
with the corresponding errors in shaded dark grey bands. 



were estimated by implementing a Monte Carlo integration 
technique. 

As mentioned in the previous paragraph, we have intro- 
duced two probability functions. One of them characterising 
halo shape distribution and the other describing the align- 
ment between dark matter haloes. Each of these functions 
has two free parameters which allow to quantify the influ- 
ence of halo alignment and halo shape on the anisotropic 
cross-correlation function. 

As a first test of our algorithm, we have compared 
the isotropic cross-correlation functions obtained from a nu- 
merical simulation with the results of the anisotropic halo 
model averaging on a spherical volume. We have found that 
taking the triaxial nature of dark matter haloes into ac- 
count improves at least %15 the predictions of the stan- 
dard halo model, as can be seen on figure [3] This improve- 
ment is more noticeable on scales corresponding to the 1- 
halo regime. T h is is in the same direction as pointed by 
ISmith fc Watts! (|2005l ) in the sense that triaxial averaging 
of halo shapes produce small but noticeably deviation from 
the classic model. Further comparison in a quantitative way 
are difficult to perform since we compute anisotropic cross- 
correlation functions whereas these authors only calculate 
the isotropic power spectrum. 

We have compared the predictions of our model with 
the results obtained from a numerical simulation, and we 
have estimated the best fitting parameters of the alignment 
and shape probability functions by means of a maximum 
likelihood method. We have found that our model is able 
to reproduce the numerical measurements over a wide range 
of scales, particularly in the 2-halo regime. Moreover, the 
model parameters obtained by fitting these numerical re- 
sults recover, as expected, the well known mass dependence 
of halo shapes and the alignment of dark matter haloes with 
the surrounding structure. The parameters C and D de- 
scribe, in a intuitive way, the effect on the cross-correlation 
function of the possible alignment configuration between 
haloes at large scales. On the other side, the parameters 
A and B reflect the halo triaxiality at small scales, never- 
theless the simplicity of the adopted profile does not allow 
a fit as good as in the case of the 2-halo regime. 

In a forthcoming paper we shall consider and develop 
the presented model in order to compare theoretical pre- 
dictions with anisotropic correlation functions measured on 
galaxy catalogues. 
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